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I show that in Bose Glass superconductor with high j c and at low T the magnetization relaxation 
(S), dominated by quantum tunneling, is oc \f] c , which crosses over to the conventional classical 
rate oc T/j c at higher T and lower j c , with the crossover T* ~ fj" 1 ■ I argue that due to interactions 
between flux lines there exist three relaxation regimes, depending on whether B < B « Bq, 
B > B$, corresponding to Strongly-pinned Bose Glass (SBG) with large j C 2, Mott Insulator (MI) 
with vanishing S, and Weakly-pinned Bose Glass (WBG) characterized by small j c i- I discuss 
the effects of interactions on j c and focus attention on the recent experiment which is consistently 
described by the theory. 
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In high T c superconductors flux line vortices play a ma- 
jor role in determining dissipation effects, and as a result 
of strong fluctuations a fairly passive Abrikosov lattice is 
replaced by a variety of new phases. At high tempera- 
tures and fields, thermal fluctuations melt the Abrikosov 
lattice over a large portion of the B-T phase diagram re- 
sulting in a strongly dissipating vortex liquid [1]. On the 
other hand, quenched disorder, at low T and B, leads to 
a variety of proposed glassy phases with Vortex Glass [2] 
corresponding to point disorder such as oxygen vacancies 
and interstitials, and Bose Gass [3] resulting from the in- 
troduction of columnar defects into the superconductor 
[4]- 

In the presence of both point and correlated disorder 
there is a variety of regimes with very different mech- 
anism for flux line motion, depending on the relative 
strength of disorder, temperature, and current, many of 
which have been lucidly discussed in Refs. [2,3]. In this 
Letter I focus on a narrower range of behavior, relevant 
to the specific magnetization-relaxation experiment that 
will be discussed below [5] . I study vortex motion at high 
currents j ss j c and very low temperatures, deep in the 
Bose Glass regime, by adapting the quantum creep the- 
ory [6] to dissipation in the presence of columnar defects. 
I show that quantum relaxation can be distinguished 
from the classical regime not only by the absence of tem- 
perature dependence but also by the counter-intuitive in- 
crease of the magnetization relaxation rate with j c , which 
I find to be oc \f] c in the quantum dissipation-dominated 
regime. 

Studies of magnetization relaxation have provided us 
with important clues into the mechanism of flux line mo- 
tion in the presence of currents and disorder. There ap- 
pears a natural separation of two time scales. An expo- 
nentially fast relaxation of B to the Bean critical state 
profile [7] takes place as long as local currents exceed the 
critical current j c , determined by the strength of pinning 
of vortices, followed by a very slow, typically logi relax- 



ation, the regime that I focus on in this work. 

Our understanding of the motion of vortices goes back 
to the Anderson-Kim (AK) model [8]. The flux lines 
are treated as independent particles thermally activated 
over finite impurity barriers U, moving down a wash- 
board potential representing disorder with the average 
slope proportional to the external current j. Because the 
barriers are finite and flux lines therefore can be ther- 
mally activated over the local barrier, this theory pre- 
dicts a finite linear flux flow resistivity for j 3> j c and 
for j <C j c (although strongly suppressed by exp (-U/T) 
factor). Near j c the I-V characteristic is highly nonlin- 
ear with V(j) ~ exp(j/j c — 1)^. Recently, however, it 
has been appreciated that for j < j c in glass phases the 
vortex motion is via vortex loop nucleation over barri- 
ers. This leads to effective barriers U(j) ~ (l/j) M that 
diverge as j — > 0, where /x is determined by the details of 
the glass phase and by the range of j. [2,3,9] Although 
these theories disagree with AK theory at low j, predict- 
ing vanishing linear resistivity, all the classical models of 
vortex motion, for j < j c , lead to the following form of 
magnetization relaxation out of the Bean state, 
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which leads to classical magnetization relaxation rate 
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An experimentally verifiable signature of this classical 
prediction is that S c vanishes linearly with T and is in- 
versely proportional to j c = Uc/<j) b, where <f> = hc/2e, b 
is the effective diameter of the columnar defect and both 
U and b are likely to be renormalized by the thermal 
fluctuations and the interactions. 

In numerous recent experiments at very low T, how- 
ever, the creep rate appears to extrapolate to a finite 
value at T = [10] which has been taken as evidence 
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of quantum tunneling of vortices out of pinning sites 
[6,11,12]. Using instanton methods it can be shown [13] 
that the thermal depinning rate exp(— U/T) is replaced 
by the tunneling rate exp(— Se/K) where Se is the Eu- 
clidean action evaluated at the stationary field configu- 
ration. In analogy to Eq.2, the quantum magnetization 
relaxation is given by [6] 



S ~ — 

DE 
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an appealing result since h analogously to T controls 
quantum fluctuations and Se plays the role of a quantum 
barrier in 3 + 1-space-time dimensional field theory. 

A transparent way to understand the form of Se is to 
look at the equation of motion for the flux line 
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where f(z,t) is the displacement of the flux line trans- 
verse to z-axis, M is the effective vortex mass, r\ and 
a are the viscous and Hall drag coefficients, v s is the 
local superfluid velocity, and F is the vortex free en- 
ergy. In the above equation, the first two terms are 
the nondissipative inertial and Hall forces, and the third 
term is the Bardeen-Stephen dissipative friction force, 
with 7] ss 4> H C 2/p n c 2 . It is not difficult to include all 
three terms in the analysis, however in high T c materials 
for j « j c and at long times I expect that the viscous 
drag term will dominate and I therefore focus on this 
term [14]. The effects of the nondissipative inertial and 
Hall terms have been studied in Refs. [6,11] for point dis- 
order and by me in the case of columnar defects. The 
inertial term leads to a T and j c -independent contribu- 
tion to S q , while the Hall term (by dimensional analysis) 
has the same scaling as the viscous contribution. The 
functional integration over r (analytically continued to 
imaginary time and transformed to frequency domain) of 
the resulting simplified equation leads to the Euclidean 
action 
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where the second term is the elastic energy, U(r) is the 
pinning energy, and the last term is the contribution due 
to the Lorentz force, //, = <\> z x j/c, acting on the flux 
line. To be specific I consider the regime in which the co- 
herence length £ « 15 A is significantly smaller than the 
diameter of the columnar defect b « 70A. In this case, 
U(f) = -U « —e a = -(</> /47tA) 2 , for r < b and van- 
ishes for a flux line outside of the columnar defect. For 
B > </> /A 2 , when the average spacing between flux lines 



a = (tpc/B) 1 / 2 3> A, the tilt modulus is well approxi- 
mated by e; ss 7e Q , with 7 = (M±/M z ) log(a/£) describ- 
ing material anisotropy. The viscous drag term (77) in 
Eq.5, is identical to the contribution to the effective ac- 
tion introduced previously by Caldeira and Leggett [15] 
to account for dissipation in a quantum system. 

To compute the quantum relaxation rate S q , I first 
determine the space-time field configuration of the tun- 
neling vortex loop, characterized by the l±, l z space 
length scales and time scale t , by balancing terms in 
the action. Balancing the pinning energy (U = b(j> j c /c) 
against the Lorentz energy, Ul z « j<f> lzl±/c, leads to 
l± ~ &(ic/i), which in a magnetization relaxation ex- 
periment with j « j c reduces to l± « b. The relax- 
ation is therefore dominated by half-loops and only in- 
volves individual columnar pins in each nucleation pro- 
cess [3]. Competition between the elastic and Lorentz 
energies, e x b 2 /l z « jc<Pohb/c, gives l z « 6(ji/j c ) 1/2 , 
where ji = eic/<j) b is proportional to the depairing cur- 
rent jo- Finally, balancing these free energy contribu- 
tions against the dissipation term, rjb 2 /t « t\b 2 /l 2 gives 
t Q « (r\jt\)l 2 z . Substituting this instanton configuration 
inside Se ~ £it Q b 2 jl z and using Eq.3 I obtain, 
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The above result obtained for individual columnar defect 
pinning is surprisingly similar to the case of collective 
quantum creep for point disorder, previously obtained by 
Blatter, et al. [6]. This Bose Glass result is suppressed by 
the factor of (&/£) 3 relative to the point disorder result, 
and the depairing current j is replaced by ji « j(£,/b)j - 
This counter- intuitive dependence of S q , in which the re- 
laxation rate increases with the pinning strength (larger 
j c ) is a clear and intriguing sign of quantum vortex mo- 
tion. 

The total magnetization relaxation rate S is a combi- 
nation of quantum and classical rates, which have very 
different j c and temperature dependences. Although a 
detailed crossover function is difficult to compute, much 
of the interesting information can be obtained from the 
asymptotic behavior of S 



S = a qn + a qd ^JJ c , for T < T* 
T 



a c — , for T > T* 

Jc 
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where a q d is defined by Eq.6, a c « c/b(j> (see Eq.2), and 
a qn is j c -independent contribution due to quantum non- 
dissipative inertial term that I expect to be small [6]. 
The crossover as a function of temperature from con- 
stant quantum to linear classical relaxation previously 
has been observed in several experiments. Balancing the 
quantum and classical relaxation rates in Eq.7 predicts 

3 /2 

the crossover temperature T* ~ j c ' . Furthermore, if a 
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consistent study of S as a function of j c can be made, I 
predict that quantum \ff c dependence dominant at high 
critical currents will crossover to classical l/j c behav- 
ior in samples with smaller j c , and this will occur at 
j* ~ T 2 / 3 . In addition, a crossover as a function of j c 
from a dissipative to a nondissipative regime can take 
place at low T and j c (see Eq.7). 

The dissipative quantum relaxation regime predicted 
by Eq.7 and a hint of crossover to classical relaxation as 
a function of j c might have already been observed in re- 
cent experiments on YBCO single crystals at milliKelvin 
temperatures. The samples were irradiated with heavy 
ions to create a forest of columnar defects with density 
characterized by = 4> /d 2 (where d is the average 
distance between the defects) , and the magnetization re- 
laxation has been measured [5] . In Figure 1 I displayed a 
smooth interpolation between the quantum and classical 
limits of Eq.7 and fit it to the S(j c ) data that I extracted 
from the experiment. At higher j c the data clearly shows 
a counter-intuitive increase in S as the pinning strength 
(j c ) is increased, which I could not explain by any classi- 
cal mechanism, in agreement with the dissipative quan- 
tum tunneling predicted by Eq.7. The flattening out at 
low j c is also consistent with the crossover to classical 
relaxation, as discussed above. 




jc 

FIG. 1. Relaxation data at T = 0.1 taken from Ref.5, plot- 
ted as a function of j c . The high (low) j c pair of points cor- 
responds to data in SBG, B <C -B</> (WBG, B ^> B^) regimes 
for B,p = 1 and 2 Tesla samples, respectively. The solid curve 
is a smooth interpolation between the quantum and classi- 
cal limits of Eq.7, fit to the experimental data. The inset 
shows the j c step-like behavior as a function of B predicted 
by the existence of SBG and WBG regimes, consistent with 
the experiment. 

Although this quantum relaxation calculation Eq.7 ap- 
pears to explain the qualitative feature of the intriguing 
increase of S with j c observed in the experiment of Ref. 



[5], an important puzzle remains. My calculation is only 
strictly correct for the regime where the pinning is dom- 
inated by individual defects. While this regime is appro- 
priate to the range of fields explored in the experiment 
of Ref. [5], it is difficult to simultaneously reconcile the 
observed increase of j c (pinning strength) with the in- 
creased density of pins (B^) in this single defect pinning 
regime. It is clear that in Eq.7 the length scale b and 
pinning energy U = b(p j c /c should be interpreted as the 
effective pinning length and pinning energy, renormal- 
ized by flux line interaction and temperature. However, 
it is not at all clear how these renormalized quantities 
are affected by the density of columnar defects, i.e. their 
dependence on B$. In the absence of a detailed theory 
of this dependence, I take the increase with B$ of the 
renormalized pinning energy U (j c ) appearing in Eq.7 
as an experimental fact (that requires further investiga- 
tion). Assuming that the change in density of defects 
has weaker affect on b than the observed increase of j c , 
Eq.7 provides a reasonable explanation for the counter- 
intuitive experimental observation of the increase of S 
with j c . 

I now turn to the effects of interactions between flux 
lines on the above picture. In examining the role of inter- 
actions on pinning in the Bose Glass phase I have found 
two quite distinct pinning mechanisms, for B < B^ and 
B > B^ to which I refer as Strongly-pinned Bose Glass 
(SBG) and Weakly-pinned Bose Glass (WBG), respec- 
tively. For B < Bfy the columnar defects outnumber 
the flux lines, and therefore all the vortices are strongly 
pinned by individual defects. The pinning energy in 
this case is of order U « e Q log(&/£) and j c i « Uc/4> b. 
As B is increased past B^, a competition between pin- 
ning energy and magnetic repulsion ensues, resulting in 
lowering of the effective pinning potential. The addi- 
tional flux lines can attempt to double-occupy the al- 
ready occupied columnar defects, thereby gaining con- 
densation energy U, but raising the interaction energy 
by Vi n t ~ 2e Q log(a/6). In these experiments a « 400A, 
b « 70 A, the magnetic repulsion dominates and the ad- 
ditional flux lines go into the interstitials in between the 
vortices localized on the columnar defects. Hence for 
B > B^ the weakly pinned interstitials flux lines are lo- 
calized by magnetic repulsion from the vortices strongly 
pinned by the columnar defects. Although the pinning 
potential seen by the interstitial vortices is not isotropic 
in the ab-plane, because the potential is ^-independent, 
the Bose Glass phase is expected to persist for B > B^, 
as was first argued in Ref. [3]. I expect this WBG phase, 
however, to have a significantly reduced critical current, 
jd relative to that of the SBG phase. 

The step- like behavior of j c , shown qualitatively in the 
inset of Fig.l, has necessarily experimentally verifiable 
consequences that the spatial linear field profile of the 
Bean state will be replaced by a two slope configuration 
corresponding to j c \ and j C 2 , with the change in the slope 
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occurring at x{B ( j ) ). The measurements of the hysteresis 
loop already find j c (B < B$) = j c2 >> j c (B > B$) = 
j c i in the strong support of the existence of distinct SBG 
and WBG pinning regimes. [5] 

I further observe that in SBG regime both B^ = 1 and 
2 Tesla samples studied in Ref. [5] exhibit a tempera- 
ture independent relaxation rate up to 4 Kelvin, giving 
additional evidence for quantum creep mechanism. On 
the other hand in the WBG regime, with significantly 
lower j c , the relaxation rate is approximately linear with 
temperature down to T* ss OAK , manifesting a classical 
relaxation mechanism. This is in qualitative agreement 
with Eq.7 which gives a quantum, T-independent rate for 
high j c and a classical, linear in T rate for low j c . 

Another striking observation is the strong suppression 
of S over a narrow range of B around B ss B^. The 
location of the dips in S suggests that the experiment 
is the first observation of the Mott Insulator phase first 
discussed in Ref. [3]. In this regime the number of flux 
lines matches the number of columnar defects, and due to 
magnetic repulsion between vortices the hopping rate is 
expected to be significantly reduced. Although the quali- 
tative picture is clear, a more thorough understanding of 
the interactions is needed to describe S near B$ in detail. 

I now examine the behavior of j c as a function of the 
columnar defect density B$. I note that j c is observed 
to approximately double as B^ is increased from 1 to 2 
Tesla. This is expected since higher density of columnar 
defects increases the average pinning strength by the ra- 
tio of attractive pinning area to total area (b/d) 2 oc B^. 
Flux line interactions lead to another source of B$ de- 
pendence in j c . Since in the absence of defects the 
flux lines would like to form an equally spaced lattice, 
in order to take advantage of the columnar defect, a 
flux line must displace an average distance d/2 from 
its otherwise preferred position. Hence the effective en- 
ergy gain from localizing on the defect is reduced by 
the increase in the repulsion with (A/a) 2 other lines to 
U eff « U - Se (d/2a) 2 (\/a) 2 « e (l - 5B 2 /AB^B cl ) 
(for B « B$, a* < A). This dependence predicts an 
additional slower increase in j c with B$ and leads to a 
critical field B* « 2^jB < j,B c i/8 above which the effective 
pinning strength vanishes. In the case a* > A, B* « B$ 
[3]. Turning the argument around, I therefore expect a 
critical B^ (B) density of defects for a realization of a true 
Bose Glass phase. 

In summary, I have discussed a theory of magnetization 
relaxation at ultra low temperatures in a Bose Glass su- 
perconductor, predicted a quantum to classical crossover 
as a function of j c as a new signature of quantum relax- 
ation in addition to the usual T-dependence signature, 
and have argued for the existence of SBG and WBG dis- 



tinct pinning regimes. I have examined recent experi- 
ments, and found that they can be consistently explained 
in the context of the developed picture. 

Note added: After this paper was submitted for publi- 
cation, one of the referees alerted me of the existence of 
the paper by V. Vinokur, Physica A Vol. 200, p. 384 
(1993), which also explores quantum vortex tunneling 
from a columnar defect. My findings are in disagreement 
with Vinokur's, who finds that the relaxation rate (in the 
regime considered here) is independent of j c , inconsistent 
with experiments of Ref. [5] . 
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